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During  this  period,  research  has  been  completed  in  a  number  of 
distinct  areas.  The  results  obtained,  as  well  as  the  goals  of  a  number 
of  related  ongoing  investigations,  are  detailed  below: 

1.  '  Estimation  in  the  NBU- Cbtsl  Class  of  Survival  Distributions  1  (Reneau, 


Samaniego) 

Suppose  the  survival  function  S  of  a  random  variable  X  is  a  member 
of  the  NBU- Ct0)  class.  By  definition,  S  must  satisfy  the  following 
condition: 


S(x+t0)  <  S(x)  S(tQ)  for  all  x  >  0. 


Let  CXi . XnJ  be  a  random  sample  from  S.  Define  the  indicator  function 

of  the  set  A  to  be 

/I  if  x  e  A 
A  X  ^  0  otherwise 

S  may  be  estimated  using  the  empirical  survival  function  Sn,  defined  by 

Sn(t)  ‘  n  ih  ' 

Sn(t)  is  simply  the  proportion  of  observations  in  the  sample  which  are 
greater  than  t.  Sn  is  a  strongly  consistent  estimator  of  S  which 
converges  to  S  at  an  optimal  rate.  However,  Sn  may  not  necessarily 
satisfy  condition  (1),  i.e  .  Sn  need  not  be  an  NBU-Ct0)  distribution.  Our 
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goal  is  to  find  an  estimator  that  itself  belongs  to  the  NBU-(t0}  class  and 
has  the  sane  asymptotic  optimality  properties  as  Sn. 

A 

In  this  paper,  the  estimator  S  of  S,  given  implicitly  by 


a  /  Sn(x) 

s(*>  -  \  .i„  r 


0<x<t, 


X  >  t. 


min  [ Sn(x) ,  S(x-t0)S(t0) ; 


is  studied.  The  following  results  are  obtained: 


Theorem  1:  Condition  (2)  uniquely  defines  an  NBU-£t0)  survival  function 
which  may  be  written  as  a  function  of  Sn  as  follows: 

(3)  S(x)  -  min  {  sfc(t0)  Sn(x-kt0)  }  , 

0sks[x/to] 

where  [u]  denotes  the  greatest  integer  less  than  or  equal  to  u. 

A 

Theorem  2:  S  is  the  largest  NBU-Ct0)  survival  function  which  is  less 

than  or  equal  to  Sn  for  all  x  >  0. 

Theorem  3 :  If  S  is  in  the  NBU-£t0)  class,  then  S  converges  uniformly  to 

S  with  probability  1.  i  Fe 

&T 

A 

Theorem  4:  The  rates  of  mean  square  and  almost  sure  convergence  of  S  to  S,3(j 


are  given  by: 


E|S(x)-S(x)|  “  0(l/7n) 


■  y - - - - — - — 

D* t  rlbuti  on/ _ 

vail  a  blllty  Code* 
jvivall  an<2/or 
nist  i  r.Doolal 


S (x) -S(x) |  -  o([L°S  ^  S]"») 


These  results  establish  the  fact  that  the  estimator  S,  given  implicitly 

/ 

in  (2)  and  explicitly  in  i !>  above,  is  a  strong  uniformly  consistent  NBU-  [ 

C t0 3  estimator  of  S  which  w.-erges  to  S  at  the  best  possible  (pointwise 


□  □ 


and  mean  square)  races.  Confidence  procedures  based  on  S  and  also 
discussed.  A  technical  report  on  this  work  will  appear  in  Summer,  1987. 


j 

2.  Estimation  of  Survival  Distributions  which  are  New  Better  Than  Used 

r - 

JWith  Respect  to  a  Distinguished  Set.  (Reneau,  Samaniego) 

)  ; 

For  any  set  A  c  [0,®),  we  define  a  survival  function  S  to  be  New 

Better  Than  Used  with  respect  to  the  Set  A  (NBU-A)  if 


S(x|t)  5  S(x|0)  for  all  x  >  0,  t  e  A 


where  S(x|t)  ■  P(X>x+t | x>t) .  The  condition  is  equivalent  to 


(4)  S(x+t)  <  S(x)  S(t)  for  all  x  >  0,  t  e  A  . 


This  class  clearly  contains  the  NBU  class  of  distributions,  and  it 
is  in  turn  contained  in  the  NBU-Ct)  class  for  all  t  S  A.  In  fact,  we  can 
write 


(5)  NBU-A  -  n  NBU-Ct)  . 

teA 

Recursive  estimation  formulas  such  as  (2)  above  are  studied  for  a  number 
of  special  cases.  When  A  is  a  finite  set,  the  resulting  estimators  are 
shown  to  be  strong  uniformly  consistent  estimators  of  the  corresponding 
survival  curves.  The  same  is  true  for  the  set  A  -  [k  tQ  :  k  J  1).  When 
the  set  A  contains  a  limit  point,  the  methods  used  on  the  preceeding  cases 
are  not,  in  general,  successful  in  producing  a  consistent  estimator.  They 
do  provide  useful  estimators  in  the  following  important  application. 

Consider  the  class  of  survival  functions  S  in  the  NBU-A  class  for 
which  there  exists  6  >  0  such  that  S(x)  -  1  for  0  <  x  <  6.  One  situation 
where  this  type  of  survival  function  occurs  in  practice  is  when  products 
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are  subjected  to  burn-in.  It  is  common  for  manufactured  items  to 
experience  a  disproportionate  number  of  failures  during  the  early  portion 
of  their  lifetimes.  Dhillon  and  Reiche  (1985)  list  inadequate  quality 
control  and  manufacturing  methods,  substandard  materials  and  workmanship, 
wrong  start-up  and  installation,  difficulties  in  assembly,  inadequate 
handling  methods  and  wrong  packaging  as  some  of  the  reasons  for  these 
early  failures.  One  method  to  deal  with  this  phenomenon  is  the  burn-in 
process  in  which  items  are  subjected  to  a  preliminary  screening,  and 
those  that  fail  before  a  predetermined  time  S  are  discarded.  The  survival 
distribution  S  of  the  items  which  survive  the  burn-in  period  will  satisfy 
the  condition  S(x)  -  1  for  x  <  6. 

The  survival  function  obtained  by  applying  left- truncation  to  a 
survival  function  S, 

f  1  if  0  <  x  <  S 

Sfi(x)  - 

l  S (x)  /  S  (5)  if  6  <  x 

is  itself  in  the  NBU-A  class  whenever  S  is.  Thus,  random  sampling  from  a 
NBU-A  distribution  under  fixed  left- truncation  such  as  burn-  in  yields  a 
survival  function  of  the  desired  type.  When  S  has  the  property  that  S(5) 
-  1  for  some  S  >  0,  the  recursive  approach  to  estimating  S  is  shown  to 
lead  to  a  strong  uniformly  consistent  estimator  with  optimal  convergence 
rates.  A  technical  report  on  this  work  is  in  preparation. 


conditions  for  the  lifetime  of  one  system  to  be  stochastically  larger  than 
that  of  another  system.  This  paper  has  been  submitted  for  publication. 


4 .  On  Failure  Rate  Estimation  When  Test  and  Operating  Environments 

Differ.  (Samaniego) 

Lindley  and  Singpurvalla  [2]  have  studied  exponentially  distributed 
failure  time  data  when  components  are  tested  in  two  different  environments 
and  each  component  failure  rate  in  the  second  environment  is  a  constant 
multiple  of  its  failure  rate  in  the  first  environment.  In  this  paper,  the 
maximum  likelihood  estimator  of  this  model's  parameter  vector  is  derived 
in  terms  of  the  unique  positive  root  of  a  specific  rational  function.  The 
estimator  of  component  failure  rates  is  compared  to  the  MLE  of  these 
failure  rates  when  no  scale -change  assumption  is  made.  This  paper  has 
been  submitted  for  publication. 

5.  Joint  Life-Time  Distributions.,  (Ghurye) 

>  l— 

The  class  of  distributions  having  the  lack-of -memory  property  was 
extended  to  incorporate  simple  aging  characteristics.  These  results  have 
now  been  published  in  Advances  in  Applied  Probability  (April  1987) . 

Since  such  a  distribution  depends  on  a  large  number  of  parameters, 
considerable  simplification  is  achieved  by  restricting  attention  to  a 
system  all  of  whose  components  are  identical;  in  such  a  situation,  the 
life-times  are  exchangeable  variables,  which  reduces  the  analytical 
complexity  to  the  minimum  possible  (with  the  exception  of  the  trivial 
case  of  mutually  independent  components).  We  are  then  dealing  with  a 
system  (such  as  "k-out  -■>:  •  which  has  several  operational  states,  and 
can  be  maintained  in  opt' re.  -li  condition  by  occasional  replacement  of 
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some  of  the  failed  components.  We  have  now  shown  that  the  solution  of 
this  latter  problem  depends  only  on  (more-or-less)  routine  computations. 

6.  Stochatic  Differential  Equations  for  Repair  Policies.-  (Ghurye) 

Optimization  problems  of  multi -component  systems  involve  large, 
though  finite,  state-spaces  (as  in  the  example  discussed  in  the  previous 
section) ;  such  problems  are  usually  soluble  by  specific  computations  in 
each  case,  but  not  in  terms  of  general  formulae  which  help  in  the 
understanding  of  the  basic  phenomena.  Therefore,  it  might  be  worthwhile 
to  replace  the  discrete  sets  by  continuous  idealizations;  this  approach 
has  proved  fruitful.  A  brief  description  follows. 

Let  X(t)  denote  the  state  of  the  system  at  time  t  (the  state  being 
measured  by  the  level  of  deterioration  or  damage),  r(t)  the  rate  of  repair 
being  done  and  i(t)  the  load  on  the  system  (such  as  pay- load  carried  or 
out-put  produced);  then  i  helps  S  to  increase  X,  whereas  r  decreases  it. 
Since  the  forces  affecting  the  system  come  from  random  processes,  the  X- 
process  satisfies  a  stochastic  differential  equation.  The  problem  of 
determining  an  optimal  repair-and- loading  policy  (r,2)  therefore  becomes, 
in  principle,  a  stochastic  control  problem.  But  whereas  most  of  the 
stochastic  control  literature  deals  with  continuous  processes, 
discontinuous  shock  processes  play  a  predominant  role  in  our  context. 
Types  of  results  which  have  been  obtained  are  illustrated  with  reference 
to  a  simple  model: 

Suppose  that  the  system  is  subject  to  damage  from  two  sources;  one  of 
these  is  a  simple  Poisson  process,  and  the  other  is  a  Brownian  motion 
with  standard  deviation  ar.<i  drift  parameter  equal  to  1.  Suppose  also 
that  a  damage -caus ing  stimulus  of  unit  value  produces  damage  of  magnitude 


iX.  Then  we  have  the  following  simple  stochastic  differential  equation: 

dX  -  [iX-r]  dt  +  a iX  dZ  +  iX  dY 


An  optimal  adaptive  stochatic  control  policy  which  maximizes  the  expected 
value  of  discounted  utility  up  to  a  finite  horizon  H  consists  of  f*  •  a 
constant  !0  and  r*  -  a(t)X(t),  where  a  is  an  exponential  function.  In  the 
limit,  as  a  goes  to  0,  i0  also  converges  to  0,  and  a  to  ®.  In  other 
words,  if  the  only  process  affecting  the  system  is  a  simple  Poisson  shock 
process,  then  no  differentiable  policy  is  admissible.  On  the  other  hand, 
if  we  admit  instantaneous  quantal  repairs  (replacements  of  components  are 
essentially  of  this  type),  then  there  exists  a  repair-and-  loading  policy 
which  prescribes  a  constant  load  and  guarantees  totally  risk-free 
operation  upto  the  horizon  H  (that  is  to  say,  a  policy  with  zero 
probability  of  a  breakdown).  Unfortunately,  this  policy  is  independent  of 
H,  and  is  identical  with  the  infinite-horizon  policy.  Consequently,  if 
the  system  is  facing  a  short-term  emergency,  during  which  it  is  necessary 
to  place  a  heavier  load  on  the  system  than  is  advisable  in  long-term 
usage,  this  condition  is  possible  only  at  the  cost  of  a  positive 
probability  of  total  breakdown  before  the  horizon  is  reached.  Further 
work  on  problems  of  this  type  is  in  progress. 

7 .  Positive  Dependence,  Upper  Sets,  and  Multidimensional  Partitions. 

(Whitaker,  A.R.  Sampson) 

Let  U  be  an  upper  set  contained  in  the  finite  discrete  lattice  L  - 

{1 . D(l)}  X  ...  X  (1,  ,D(p)).  Representations  for  U  are  obtained 

and  shown  to  correspond  to  certain  multidimensional  partitions  of 
integers.  It  is  shown  th.r  tor  p-3,  the  number  of  possible  upper  sets  in 


Various  other  representation  and  enumeration 


L  1S 

results  are  obtained  for  related  settings.  Considered  are  a  variety  of 
applications  to  multivariate  positive  dependence  and  various  notions  in 
statistics  and  reliability  theory.  A  paper  based  on  this  work  will 
appear  in  Mathematics  of  Operations  Research. 

8.  Estimation  of  Multivariate  Distributions  Under  Stochastic  Ordering. 
(Whitaker,  A.R.  Sampson) 

Let  F  and  G  be  the  cdf's  of  two  p-dimensional  multivariate 
distributions,  such  that  F  is  stochastically  larger  than  G.  A 
straightforward  derivation  is  given  of  the  generalized  maximum  likelihood 
estimators  of  F  and  G  based  on  the  random  samples  from  each  population. 
An  algorithmic  approach  to  computing  these  estimators  is  described  and 
motivating  numerical  examples  are  discussed.  The  special  case  when  F  and 
G  correspond  to  multivariate  ordinal  contingency  tables  is  also  presented. 
The  relationship  of  these  results  to  those  of  Robertson  and  Wright  (Ann. 
Statist.  2  (1974)  528-534)  is  considered.  A  paper  on  this  work  has  been 
submitted  for  publication. 

9 .  ^Computational  Aspects  of  Association  for  Bivariate  Discrete 

Distributions .  (Whitaker,  A.R.  Sampson) 

For  a  bivariate  discrete  probability  distribution,  P,  on  the  RxC 
lattice,  various  aspects  for  checking  association  (Esary,  Proschan  and 
Walkup  (1967))  are  considered.  A  new  algorithm  is  given  for  verifying 
whether  or  not  P  is  associated.  The  efficiency  of  this  algorithm  is 
obtained  and  compared  to  the  etficiency  of  a  simple  algorithm  based  on  the 
definition  of  association  -lien  R-C-5  for  example,  the  new  algorithm 
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requires  less  than  3%  of  the  computations  required  by  the  simple 
algorithm.  In  obtaining  these  results,  a  new  set  function  Q  is  defined  on 
all  upper  sets  in  the  lattice;  the  relationship  between  P  and  Q  sheds 
light  on  the  relationship  between  association  and  stochastic  ordering. 
Also  defined  is  a  computationally  important  subset  of  the  extreme  points 
of  an  upper  set.  A  technical  report  of  this  work  is  in  preparation. 


In  addition  to  the  work  described  above,  the  investigators  have 
revised  a  number  of  grant  supported  research  papers  in  various  stages  in 
the  publication  process.  These  include  the  following: 


"Estimating  a  Survival  Curve  Based  on  Nomination  Sampling,"  Journal  of 

the  American  Statistical  Association,  81(1986),  1034-45  (Boyles  and 
Samaniego) . 

"On  Estimating  Component  Reliability  for  Systems  with  Random  Redudancy 
Levels,"  IEEE-Transactions  on  Reliability,  to  appear  (Boyles  and 
Samaniego) . 

"On  Estimating  Population  Characterisitcs  From  Record-Breaking 

Observations,  II.  Nonparametric  Results,"  Naval  Research  Logistics, 
to  appear  (Samaniego  and  Whitaker). 

"On  Estimating  the  Reliability  of  Systems  Subject  to  Imperfect  Repair," 

Journal  of  the  American  Statistical  Association,  to  appear  (Whitaker 
and  Samaniego) . 

"Estimating  a  Survival  Curve  when  New  is  Better  Than  Used  in  Expectation, 
submitted  for  publication  (Whitaker  and  Samaniego) . 

"New  Moment  Identities  Based  on  Integrated  Survival,"  submitted  for 
publication  (Arnold,  Reneau  and  Samaniego). 
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